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Abstract 

Quasiparticle dynamics in relativistic plasmas associated with hot, weakly-coupled gauge theo- 
ries (such as QCD at asymptotically high temperature T) can be described by an effective kinetic 
theory, valid on sufficiently large time and distance scales. The appropriate Boltzmann equations 
depend on effective scattering rates for various types of collisions that can occur in the plasma. 
The resulting effective kinetic theory may be used to evaluate observables which are dominantly 
sensitive to the dynamics of typical ultrarelativistic excitations. This includes transport coeffi- 
cients (viscosities and diffusion constants) and energy loss rates. In this paper, we show how to 
formulate effective Boltzmann equations which will be adequate to compute such observables to 
leading order in the running coupling g{T) of high-temperature gauge theories [and all orders in 
1/ log g{T)~^]. As previously proposed in the literature, a leading-order treatment requires includ- 
ing both 2^2 particle scattering processes as well as effective "1 2" collinear splitting processes 
in the Boltzmann equations. The latter account for nearly collinear bremsstrahlung and pair pro- 
duction/annihilation processes which take place in the presence of fluctuations in the background 
gauge field. Our effective kinetic theory is applicable not only to near-equilibrium systems (rel- 
evant for the calculation of transport coefficients), but also to highly non-equilibrium situations, 
provided some simple conditions on distribution functions are satisfied. 
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I. INTRODUCTION 



In a hot, weakly coupled gauge theory, such as QCD at asymptotically high temperature 
where the running coupling g{T) is small, one might hope to achieve solid theoretical control 
over the dynamics of the theory. To date, however, very little has been derived about the 
dynamics of such theories at even leading order in the coupling — that is, neglecting all 
relative corrections to the leading weak-coupling behavior which are suppressed by powers 
of g. For example, hydrodynamic transport properties such as shear viscosity, electrical 
conductivity, and flavor diffusion are not known at leading order; they have only been 
calculated in a "leading-log" approximation, which has relative errors of order 1/ log(5'~^).^ 

To study transport or equilibration processes in a hot plasma quantitatively, the most 
efficient approach is first to construct an effective kinetic theory which reproduces, to the 
required level of precision, the relevant dynamics of the underlying quantum field theory, 
and then apply this kinetic theory to the processes of interest. Specifically, we would like 
to formulate an appropriate set of Boltzmann equations which will, on sufficiently long time 
and distance scales, correctly describe the dynamics of typical ultrarelativistic excitations 
{i.e., quarks and gluons) with sufficient accuracy to permit a correct leading-order evaluation 
of observables such as transport coefficients. Schematically, these Boltzmann equations will 
have the usual form. 



where / = f{x,p,t) represents the phase space density of (quasi-)particles at time t, v 
is the velocity associated with an excitation of momentum p, and C[f] is a spatially-local 
collision term that represents the rate at which particles get scattered out of the momentum 
state p minus the rate at which they get scattered into this state. The challenge is to 
understand exactly what processes need to be included in the collision operator C[f], and 
how to package them, so that the Boltzmann equation correctly reproduces the desired 
physics at the required level of accuracy. 

To compute transport coefficients or asymptotic equilibration rates, one does not actually 
need a general non-equilibrium (and non-linear) kinetic theory; it is sufficient merely to have 
Boltzmann equations linearized in small deviations away from an equilibrium state of given 
temperature T.^ But more generally one would like to formulate a fully non-equilibrium 
kinetic theory which would also be applicable to systems (such as intermediate stages of a 
heavy ion collision [3]) in which deviations from equilibrium are substantial and quantities 
such as temperature are not unambiguously defined. This will be our goal. It should be 
emphasized that all of our analysis assumes the theory is weakly coupled on the scales of 
interest. Hence, the domain of validity includes QCD at asymptotically high temperatures, 
or intermediate stages of collisions between arbitrarily large nuclei at asymptotically high 
energies [3]. It is an open question to what extent this weak coupling analysis is applicable 
to real heavy ion collisions at accessible energies. However, we believe that understanding 



^ Ref. [1] attempted to calculate the shear viscosity to next-to-leading logarithmic order [including relative 
corrections of order l/log(g~^) but neglecting l/log^(g~^) effects] but, among other things, missed the 
"1 2" collinear processes described in this paper, as well as pair annihilation/creation processes which 
contribute even at leading-log order. The latter have been previously discussed in Ref. [2] . 

^ See, for example, the discussion in Ref. [2] . 




(1.1) 
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FIG. 1: Scattering of hard particles by t-channel gauge boson exchange. The straight hnes could 
represent any type of hard charged particle, including fermions, gauge bosons, or scalars. 

the dynamics in weakly coupled asymptotic regimes is a necessary and useful prerequisite 
to understanding dynamics in more strongly coupled regimes. 

The domain of applicability of any kinetic theory depends on the time scales of the un- 
derlying scattering processes which are approximated as instantaneous transitions in the 
collision term of the Boltzmann equation. In the remainder of this introduction, we review 
the relevant scattering processes and associated time scales, describe the assumptions under- 
lying our effective kinetic theory for near-equilibrium systems in which the temperature is 
at least locally well-defined, and then discuss how the required conditions can be generalized 
to a much wider class of non-equilibrium systems. 



A. Relevant scattering processes 

Consider a QCD plasma at sufficiently high temperature T so that the effective couphng 
g(T) is small. ^ Quarks and gluons are well-defined quasiparticles of this system. The typical 
momentum (or energy) of a quark or gluon is of order T; this will be referred to as "hard." 
The number density of either type of excitation is O(T^), so that the total energy density is 
O(T^). Hard quarks and gluons propagate as nearly free, nearly massless excitations. Their 
dispersion relations receive thermal corrections which look just like an effective mass,^ 

e(p) = vW<, (1-2) 

with the thermal masses for both quarks and gluons being 0{gT) in size. Hence, all but 
a parametrically small 0{g'^) fraction of excitations [those with 0{gT) momenta] are ultra- 
relativistic and travel at essentially the speed of light. If p denotes the unit vector in the 
direction of p then, for hard excitations, 

v{p)^^^p[l-0{g')]. (1.3) 

Scattering processes in the plasma will cause any excitation to have a finite lifetime. 
Imagine focusing attention on some particular excitation with a hard 0{T) momentum. 



^ For convenience of presentation, we will use the language of QCD throughout this paper. Despite this, 
all our discussion and results apply equally well to high temperature electroweak theory and, except for 

a few comments about non-pcrturbativc non-AbeUan g^T scale physics, also to high temperature QED. 
^ This simple form holds, up to yet higher-order corrections, provided the momenta of the excitation is 
large compared to the thermal mass, |p| 0{gT). Our rrith is often called the asymptotic thermal mass. 
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TABLE I: Parametric dependence of various length scales for a weakly-coupled ultrarelativistic 
equilibrium plasma. Estimates for mean free paths apply to typical (hard) particles, with 9 and q 
denoting the deflection angle and momentum transfer, respectively. The non-perturbative magnetic 
physics scale of {g^T)'^ for colored fluctuations applies only to non-Abelian gauge theories. 



What determines the fate of this quasiparticle? One relevant process is ordinary Coulomb 
scattering, depicted in Fig. 1. If the particle of interest scatters off some other excitation in 
the plasma with a momentum transfer g, then the direction of the particle can change by 
an angle 9 which is 0{\q\/T). The differential scattering rate is 

d^-9'T^^j-9'T^i,. (1-4) 

This form holds provided q ^ 0{gT). Below this scale, Debye screening (and Landau damp- 
ing) in the plasma soften the small angle divergence of the bare Coulomb interaction. 

Consequently, the rate for a single large angle scattering with 0{T) momentum transfer is 
0{g'^T), while the rate for small angle scattering with 0{gT) momentum transfer is 0{g^T). 
For later use, let Tg = l/{g^T) denote the characteristic small angle mean free time, and 

= l/{g'^T) the characteristic large angle mean free time,^ also known as the transport 



^ Actually, because small angle scatterings are individually more probable than large angle scatterings, a 
particle can undergo Ng ^ small angle scatterings, each with 6 = 0{g), during a time of order 

NgTg ~ l/{g'^T) — the same as the time for a single q = 0(T) scattering. A succession of this many 

1 /2 

uncorrelated small angle scatterings will result in a net deflection of order Ng 9 — 0{1). Hence, a large 
deflection in the direction of a particle is equally likely to be the result of many small angle scatterings 
or a single large angle scattering. In fact, the multiplicity of possible combinations of scatterings with 
angles between 6 ^ g and ^ ^ 1 leads to a logarithmic enhancement in the large angle scattering rate, 
or a logarithmic decrease of the large angle mean free path, so that ^ [g^nog(l/9)]^^. In this paper, 
we will ignore such logarithmic factors when making parametric estimates and simply write ^ 
Nevertheless, it is important to keep in mind that contributions to r* come from the entire range of 
scatterings from 6 g to 9 1, which corresponds to momentum transfers from q gT to g ^ T. 
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FIG. 2: t-channel diagram for qq ^ gg. 

mean free time. Neither time has a precise, quantitative definition; these quantities will only 
be used in parametric estimates. The small angle mean free time Tg is [to within a factor of 
0{logg~^)] the same as the color coherence time of an excitation — this is the longest time 
scale over which it makes sense to think of an excitation has having a definite (non-Abelian) 
color [4-7]. These scales arc summarized in Table I. 

It is also important to consider processes which change the type of an excitation. Consider, 
for example, the conversion of a quark of momentum p into a gluon of nearly the same 
momentum by the soft qq gg process, depicted in Fig. 2, with momentum transfer q ~ gT. 
The mean free path for this process (or its time-reverse) is 0[l/{g^T)], just hke the large 
angle scattering time r*. As a result, such t-channel quark exchange processes are equally 
important as gluon exchange for quasiparticle dynamics, and must be correctly included 
even in leading-log evaluations of transport coefficients [2]. 

Crossed s-channel versions of Figs. 1 and 2, namely quark-antiquark annihilation and 
creation via a single virtual gluon, and gluo-Compton scattering, also proceed at 0{g'^T) 
rates. Consequently, these processes must also be included in a leading-order treatment of 
quasiparticle dynamics.^ Henceforth, whenever we refer to 2 ^ 2 particle processes, we will 
mean all possible crossings of Figs. 1 and 2 in which two excitations turn into two excitations. 

In addition to the 2 2 particle processes of Figs. 1 and 2, hard quasiparticles in 
the plasma can also undergo processes in which they effectively split into two different, 
nearly coUinear, hard particles. Such processes cannot occur (due to energy-momentum 
conservation) in vacuum, but they become kinematically allowed when combined with a 
soft exchange involving some other excitation in the plasma. A specific example is the 
bremsstrahlung process depicted in the upper part of Fig. 3. Here, one hard quark undergoes 
a soft (q ~ gT) collision with another and splits into a hard quark plus a hard gluon, each 
of which carry an 0(1) fraction of the hard momentum of the original quark and both of 
which travel in the same direction as the original quark to within an angle of 0{g). The 
mean free path for this process, as well as the near-coUinear pair production process also 



Yet softer momentum transfers are screened in the plasma and do not affect r* at the order of interest. 
For instance, the mean free path for very-small angle scattering with 6 (corresponding to momentum 
transfers of order fif^T) is only Tg2 ~ l/g^T and is not enhanced over Tg. Over the time r, ^ 1/g^T, there 

can thus be only Ng2 ~ 1/g'^ independent very-soft scatterings, which will only contribute a net deflection 
of order A9 ^ (iVg2)^/^g^ ^ g to the 0(1) deflection caused by other processes. Consequently, the large 
angle scattering time r* is insensitive to non-perturbative magnetic physics in the plasma associated with 
very soft momentum transfers of order g'^T. 

These s-channel processes do not have the logarithmic enhancement mentioned in footnote 5, and so do 
not contribute to transport coeflacients at leading-log order, but do contribute at next-to-leading log order. 
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FIG. 3: Simplest examples of near-collinear processes which contribute at leading order. The upper 
diagrams depict hard, nearly collinear gluon bremsstrahlung accompanying a soft gluon exchange 
between hard particles. The lower diagrams show the conversion of a hard gluon into a nearly 
collinear quark-antiquark pair when accompanied by a soft exchange with another excitation in 
the plasma. 

shown in Fig. 3, turns out to be 0[l/{g^T)], which is once again the same order as the large 
angle scattering time r*.'' 

There is an important difference between these near-collinear processes and the 2^2 
particle processes of Figs. 1 and 2. The intermediate propagator in the diagrams of Fig. 3 
which connects the "splitting" vertex with the soft exchange has a small 0{g^T'^) virtuality. 
Physically, this means that the time duration of the near-collinear processes (also known as 
the scattering time or formation time) is of order 1/g^T, which is the same as the mean 
free time Tg for small-angle elastic collisions in the plasma. As a result, additional soft 
collisions are likely to occur during the splitting process and can disrupt the coherence 
between the nearly collinear excitations. This is known as the Landau-Pomeranchuk-Migdal 
(LPM) effect. Since the time between soft scatterings is comparable to the time duration 
of the emission process, multiple soft scatterings cannot be treated as independent classical 
events, but must be evaluated fully quantum mechanically. In other words, including the 
interference between different N + 1 ^ N + 2 amplitudes, such as depicted in Fig. 4, is 
required to evaluate correctly the rate for near-coUinear sphtting at leading order. We will 
refer to these processes collectively as "1 — > 2" processes where the 1 — > 2 refers to the 
nearly collinear splitting particles and the quotes are a reminder that there are other hard 
particles participating in the process via multiple soft gluon exchanges. Of course, the inverse 
"2 — > 1" processes, in which two nearly collinear excitations fuse into one, are also required 
for detailed balance. The evaluation of the rate of these "1 <-> 2" near-coUinear processes 



See, for example, Ref. [8] as well as the discussion of the closely related case of photon emission in 
Refs. [9, 10]. More careful analysis shows that the rates of these near-collinear processes do not have the 
logarithmic enhancement of the large angle scattering rate discussed in footnote 5. Therefore these near- 
collinear processes do not contribute to transport coefficients at leading-log order but must be included 
at next- to-leading log order. 



6 




FIG. 4: Examples of interference terms involving multiple scatterings that also contribute to hard- 
gluon bremsstralilung at leading order. 

in an equilibrium ultrarelativistic plasma, complete to leading order, is discussed in Ref. 
[8],^ which derives a two-dimensional linear integral equation whose solution determines the 
leading order rate. 

Having realized that near-collinear "1 ^ 2" processes are just as important to the fate 
of a quasiparticle as 2 2 particle processes, one might wonder if there are any further 
relevant processes which occur at an 0{g^T) rate and hence compete with the processes just 
discussed. We will argue in section VI B that this is not the case. 

In summary, a typical hard excitation travels a distance of order r^, ~ 1/ ((?^T) before it 
either experiences a large angle scattering, converts into another type of excitation, splits 
into two nearly coUinear hard excitations, or merges with another nearly collinear excitation. 
These different "fates" are all comparably likely [up to factors of 0{logg^^)]. During the 
transport mean free time r*, the excitation will experience many soft scatterings each of 
which can completely reorient the color of the quasiparticle, but only change its momentum 
by a small 0{gT) amount. 



B. Kinetic theory domain of validity — near-equilibrium systems 

A near-equilibrium system is one in which the phase space distribution can be written as 
an equilibrium distribution n{p) plus a small perturbation 6f{x,p). However, we wish to 
include situations where the system is close to local rather than global equilibrium, so that 
equilibrium parameters like temperature may vary slowly with x. We will therefore define a 
near-equilibrium system to be one where distribution functions can be written in the form 

f{x,p) ^ n{p;T{x),u{x), ii{x)) + 5f{x,p), (1.5) 

where n{p; T, u, /i) denotes an equilibrium distribution (Bose or Fermi, as appropriate) with 
temperature T, flow velocity u, and optionally one or more chemical potentials /x, provided: 
(i) the parameters T{x), u{x), and iJ,{x) do not vary signiflcantly over distances (or times) 
of order 1/g'^T, the large angle mean free path,^ and (ii) throughout the system, 5f -C /.^° 
In the case of near-equilibrium systems, we can now summarize the conditions upon which 
our effective kinetic theory is based. 



See also Refs. [9, 11] and references therein. 

This condition should be understood as applying in a frame which is slowly moving relative to the local 

rest frame of the fluid at point x. 

More precisely, Sf must make a negligible contribution to the integrals (1.6) and (1.7) below. This ensures 
that the local-equilibrium part of the distribution dominates screening and soft scattering processes. 
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1. As already indicated, we assume that the theory is weakly coupled on the scale of 
the temperature, g{T) -C 1. Consequently, there is a parametrically large separation 
between the different scales shown in Table I. 

2. We assume that all zero-temperature mass scales {i.e., Aqcd and current quark masses) 
are neghgible compared to the 0{gT) scale of thermal masses. 

3. Any kinetic theory can only be valid on time scales which are large compared to the 
duration of the scattering processes which are approximated as instantaneous inside the 
collision term of the Boltzmann equation. Since the formation time of near-collinear 
splitting processes is order 1 / (g'^T) , this means we must assume that the space-time 
variation of the deviation from local equilibrium 6f{x,p) is small on the scale of 
l/{g'^T). [Note that we've already assumed a stronger condition on the spacetime 
variation of the local equilibrium part of distributions, n{p;T{x),u{x), ijl{x)).] 

4. We will assume that hard particle distribution functions Sf{x,p) have smooth de- 
pendence on momentum p and do not vary significantly with changes in momentum 
comparable to the 0{gT) size of thermal masses. [Again, this is automatic for the 
local equilibrium part of distributions.] This condition will allow us to simplify our 
treatment of distribution functions for near-collinear "1 2" processes. 

5. Finally, we assume that the observables of ultimate interest are dominantly sensitive 
to the dynamics of hard excitations, are smooth functions of the momenta of these ex- 
citations, do not depend on the spin of excitations, and are gauge invariant. Excluding 
spin-dependent observables is primarily a matter of convenience, and will allow us to 
use spin-averaged effective scattering rates. Considering color-dependent observables 
would be senseless in an effective theory which is only applicable on spatial scales large 
compared to the 1 / (g^T) color coherence length. 

C. Kinetic theory domain of validity — non-equilibrium case 

It is impossible for any effective kinetic theory to be valid for all possible choices of phase 
space distribution functions. There will always exist sufficiently pathological initial states 
which are outside the domain of applicability of any kinetic theory. Describing the domain 
of applicability therefore requires a suitable characterization of acceptable distribution func- 
tions. In a general non-equilibrium setting, we will define "reasonable" distribution functions 
to be those which support a separation of scales similar to the weakly-coupled equilibrium 



For hot electroweak theory, the required condition is that the Higgs condensate be small compared to the 
temperature, v{T) <^ T, so that the condensate-induced mass of the VF-boson is negligible compared to 

its thermal mass. 

If this condition is not met, distribution functions become density matrices in spin space, and the Boltz- 
mann equation must be replaced by a generalization known as a Waldmann- Snider equation [12]. (See 
also the discussion in Ref. [4] for the analogous case of color-dependent density matrices.) We assume 
spin-independence to avoid needlessly complicating this paper, and because we are not aware of physically 
interesting problems in hot gauge theories where strong spin polarization is relevant. 
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case. That is, the momenta of relevant excitations must be parametrically large compared 
to medium- dependent corrections to dispersion relations or to the inverse Debye screening 
length. Furthermore, the phase space density of excitations must not be so large as to drive 
the dynamics, on these scales, into the non-perturbative regime. 

There are potentially many more relevant scales one may need to distinguish in a non- 
equilibrium setting. In particular, in place of the single hard scale T that characterizes 
the momenta of typical excitations in equilibrium, one may need to consider at least three 
relevant scales, which we shall refer to as the typical momenta of (i) primaries, (ii) screeners, 
and (iii) scatterers. By "primaries," we mean the particles whose evolution we are explicitly 
interested in following with the Boltzmann equation. These might be the particles which 
dominate the energy of the system (as in the intermediate stages of bottom- up thermalization 
[3]), they might be particles which dominate transport of some charge of the system we 
are following, or whatever. These excitations may or may not overlap with the next two 
categories. The "screeners" are those particles whose response to electric and magnetic 
fields dominates the screening effects and hard thermal masses in the medium. That is, 
they have momenta comparable to the momentum scale at which hard thermal-loop (HTL) 
self-energies receive their dominant contribution. (Explicit formulas for HTL self-energies 
will be reviewed in Section IV.) Finally, the "scatterers" are those particles which generate 
the soft background of gauge fields off of which the primaries scatter in soft collisions. As a 
pictorial example, consider Figs. 1, 3, and 4. The lines entering on the top are primaries, the 
lines entering on the bottom are scatterers, and the screeners are not shown explicitly but 
are the particles in the hard thermal loops that are implicitly summed in the soft, exchanged 
gluon lines. 

To formulate the conditions under which our Boltzmann equations will be valid, it will be 
helpful to define the following two integrals involving the distribution function of a particular 
species (gluon, quark, or antiquark) of quasiparticles. 



As usual upper signs refer to bosons and lower signs to fermions. In equilibrium, the ratio 
X/ J is precisely the temperature T. 

Corrections to the ultrarelativistic dispersion relation for quasiparticles will be seen to 
involve the quantity J{x)] the square of medium-dependent effective masses will be of order 
g"^ J . So the momentum dominating this integral defines the typical momentum of screeners 
which we will denote by Pscreen, and 



is the characteristic size of effective masses. This reduces to the gT scale in the near- 
equilibrium case. Momenta large compared to mcff will be referred to as "hard." 

The integral X(x) will appear as an effective density of scatterers which generate soft 
background gauge field fiuctuations. The momentum scale which dominates this integral 
is, by definition, the characteristic momentum of scatterers, Pscatter- The mean free time 




(1.6) 




(1.7) 
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between small-angle scatterings (with momentum transfer of order mefr) scales as 
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rsoft-^. (1.9) 

This generalizes the l/{g^T) small angle mean free time in the near-equilibrium case. 

Let Pprimary dcnote the characteristic primary momenta of interest, and let Phard denote the 
minimum of Pscatter, Pscreen, and Pprimary As we shall rcview in section VB, the formation time 
of near-coUinear bremsstrahlung or annihilation processes involving a primary excitation 
behaves (up to logarithms) as 

^form - J", , (1.10) 

where 

Pprimary 



-^form 



1 + 



1/2 

(1.11) 



represents the typical number of soft collisions occurring during a single near-collinear 
"1 ^ 2" process. (In other words, A^fomi is parametrically either 1 or [ppnmaxy/ {n^eS '^soh)Y^'^ : 
depending on whether ^primary is small or large compared to m^g Tsoft-) 

Since distribution functions may vary in space or time, all these scales are really local 
and may vary from point to point in spacetime. However, we will assume that distribution 
functions are suitably slowly varying in space and time. The following discussion should be 
understood as applying in the vicinity of any particular point x in the system. 

The key assumptions we will make are: 

1. The species dependence of the above scales is not parametrically large. In particular, 
ratios of the effective masses of different species are 0{1). This assumption is a matter 
of convenience, but relaxing it would significantly complicate the discussion. 

2. The momenta of primaries, scatterers, and screeners are all large compared to medium- 
dependent effective masses. A large separation of scales is essential to our analysis, 
and implies that all relevant excitations are highly relativistic. To make parametric 
estimates, we will assume that 

meff/phard < 0(^") (1.12) 

for some positive exponent a. [For instance, meff/phard was 0{g) in our previous 
discussion of near-equilibrium physics.] 

3. The effective mass rrics must be large compared to the small-angle scattering rate t^"^, 
as well as to zero-temperature mass scales (Aqcd and quark masses). 



Initial and final state factors for the scatterer appear in this expression inside the integral X, given by 
Eq. (1.7). Some readers may wonder whether there should be analogous initial and final state factors for 
the primary excitation. To be more precise, Tgoft represents the inverse of the contribution to the thermal 

width of quasiparticles due to soft scattering. This part of the width docs not contain statistical factors 
for the primaries. The scale Tgoft characterizes the time scale for the decay of quasiparticle excitations. 
Equivalently, it is the relaxation time of fluctuations in the occupancy of a hard mode due to soft scattering. 
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4. All distribution functions have negligible variation over spacetime regions whose size 
equals the formation time tform of near-coUinear processes involving primary excita- 
tions. This assumption is essential for our treatment of near-coUinear processes. 

5. Distribution functions of scatterers and screeners do not vary significantly with para- 
metrically small changes in the direction of propagation of excitations. For example, 
distributions cannot be so highly anisotropic that the directions of screeners or scatter- 
ers all lie within an 0{g) angle of each other. This simplifies the analysis by allowing 
us to ignore angular dependence when making parametric estimates. We will actually 
find that a stronger limit on the anisotropy of screeners is needed in order to prevent 
the appearance of soft gauge field instabilities [13] whose growth can lead to violations 
of the preceding spatial smoothness condition. Discussion of such instabilities will be 
postponed to section VI C. 

6. Distribution functions, for hard momenta, have smooth dependence on momentum 
and do not vary significantly with O(mefj) changes in momentum. [This assumption 
is implicit in various non-equilibrium HTL results which we will take from the liter- 
ature.] For momenta of order Pprimary; we further assume that distributions do not 

1 /2 

vary significantly with an 0{Nf^j.^mes) change in momentum, which represents the 
typical total momentum transfer due to the A^form soft collisions occurring during a 
near-coUinear "1 <-> 2" process. 

7. Distribution functions are not non-perturbatively large for momenta p > O(meff). 
Specifically, for bosonic species we will assume that f{p, x) is parametrically small 
compared to 1/g'^. (Fermionic distributions can never be parametrically large.) As 
discussed below, this inequality is actually a consequence of condition 2. 

8. Distribution functions are not spin-polarized. Once again, this condition is a matter of 
convenience. But the consistency of this assumption is now a non-trivial issue since we 
are not requiring distribution functions to be isotropic in momentum space. In a gen- 
eral (rotationally invariant) theory, it is quite possible for a medium with anisotropic 
distribution functions to generate medium-dependent self-energy corrections which 
lead to spin or polarization dependent dispersion relations, implying spin-dependent 
propagation of quasiparticles. In such a theory, initially unpolarized distribution func- 
tions would not remain unpolarized. This point will be discussed further in section 
IV, where we will see that in hot gauge theories (at the level of precision relevant for 
our leading-order kinetic theory) unpolarized but anisotropic distribution functions do 
not generate birefringent quasiparticle dispersion relations. 

9. Distribution functions are color singlets. Attempting to incorporate colored distribu- 
tion functions would be inconsistent, since the color coherence time is comparable or 
shorter than the formation time of the "1 2" processes which will be treated as 
instantaneous in this effective kinetic theory. 

10. Observables of interest are dominantly sensitive to the dynamics of excitations with 
momenta of order Pprimary, are smooth functions on phase space, are gauge invariant, 
and are spin independent. 



11 



The most important restrictions are the scale-separation condition 2 and the effective mass 
condition 3. They have numerous consequences, including condition 7, as explained below. 

Let fp denote the average phase space density for excitations with momenta of order p. 
First note that the definition of Pscreen as the momentum which dominates the integral 
(1.6) implies that 

^screen 

for any momentum p. For momenta p which are parametrically difi^erent from Pscreen, the 
above inequality ( ^ ) is actually a strong inequality (<IC) (or else Pscreen will not be the scale 
which dominates J). From (1.13) and definition (1.8), one has 

<my/. (1.14) 

For hard momenta p > Phard, condition 2 [mlg/p'^ < 0{g'^°')] thus imphes that 

fp < 0{g'--'), (1.15) 

which is a strengthened version of condition 7. For soft momenta, such as p ~ ttLcs, which 
are parametrically smaller thanpscatter by condition 2, the inequality (1.14) becomes strong, 
p^ fp <^ mlff/g'^. We can then generally conclude that for momenta p > iries, the phase 
space density is perturbative, 

f, « 0{g-') . (1.16) 

In order for any Boltzmann equation to be valid, the duration of scattering events (which 
are treated as instantaneous in the Boltzmann equation) must be small compared to the typ- 
ical time between scatterings. For 2^2 collisions, the largest relevant scattering duration 
is l/rries [for soft scatterings with 0{mes) momentum transfer] and the smallest relevant 
mean free time is the small-angle scattering time (1.9). Hence, the condition 3 requirement 
that meflf Tsoft ^ 1 is needed for the validity of kinetic theory. If the density of scatterers is 
not parametrically small, so that /p^^attcrfl ^ ipscattor] is 0(1) or larger, then this inequality 
automatically holds since Eqs. (1.13) and (1.15) imply^'^ 



meff Tsoft ~ ~ 5 Z {g /p,catteJ ' ^0{g ) . (1.17) 



7^3/2 ^3/2 3 
rvj — 

9^ gfl scatter ^scatter 

But the more general condition 3 determines the limit of applicability of kinetic theory for 
dilute systems. 

In order to regard excitations as having well-defined energies and momenta, their de 
Broglie wavelengths must also be small compared to the mean time between scatterings. The 
longest relevant de Broglie wavelength is l/phard, so the validity of kinetic theory requires 
that Phard Tsoft ^ 1- But this automatically follows from conditions 2 and 3. 

Condition 6, requiring smoothness of distribution functions in momentum space, prevents 
applications of our effective theory to cold degenerate quark matter. This condition implies 
that the temperature (for near-equilibrium systems) must be large compared to gpiermi and 



This estimate assumes that some of the relevant scatterers are bosons, which might have parametrically 
large distributions. If the only relevant scatterers are fermions, then one finds meflf Tsoft ^ 0{l/g). 
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hence, for weak coupling, lies far outside the temperature region in which color supercon- 
ductivity occurs. 

The remainder of this paper is organized as follows. Section 11 presents the structure 
of the effective kinetic theory. The effective scattering amplitudes characterizing 2^2 
processes are discussed in section 111. These quasiparticle scattering amplitudes depend 
on medium-dependent self-energies, which are the subject of section IV. The appropriate 
formulation of effective transition rates for near-collinear "1 ^ 2" processes is described in 
section V. Section VI discusses the validity of our effective kinetic theory at greater length, 
including possible double counting problems in our effective collision terms, and potential 
contributions of omitted scattering processes. We argue that neither of these concerns are an 
issue. This section also examines the possible appearance of instabilities in soft (p ~ mefr) 
modes of the gauge field, and briefly mentions open problems associated with extending our 
effective kinetic theory beyond leading order. Two short appendices follow. One summarizes 
simplifications to the formulas in the main text that can be made in the case of isotropic 
distribution functions, and the other discusses the connection between the formulas for 1 <-> 2 
scattering presented in section V and the results for the total gluon emission rate discussed 
in Ref. [8]. 



II. THE EFFECTIVE KINETIC THEORY 

Our effective kinetic theory will include all 2 ^ 2 processes as well as effective coUinear 
"1 2" processes. The Boltzmann equations are, 

{dt-rp- V.) /.(x, p, t) = - C;^-^" [/], (2.1) 

where the label s denotes the species of excitation (gluon, up-quark, up-antiquark, down- 
quark, down-antiquark, etc.). Since we have assumed that distributions are not spin or color 
polarized, we do not decorate distribution functions with any spin or color label. However it 
should be understood that fs{x,p,t) represents the phase space density of a single helicity 
and color state of type s quasiparticles.^^ 

Schematically, the overall structure of our Boltzmann equations is similar to that outlined 
by Baier, Mueller, Schiff, and Son [3] in their treatment of the late stages of their "bottom- 
up" picture of thermalization in heavy ion collisions, but our formulation of the details of 
the collision terms will be guided by our goal of providing a treatment which is complete 
at leading order. As will be discussed below, this requires a consistent treatment of both 
screening and LPM suppression of near-collinear processes. 



Since our effective theory describes typical hard excitations, gluons are to be regarded has having only two 
transverse polarizations. There is a longitudinal branch of the gluon dispersion relation in a hot plasma, 
but the spectral density of longitudinal excitations is exponentially small for hard momenta. Hence, these 
collective excitations may be completely ignored for our present purposes. 
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The elastic 2 <-> 2 collision term for a given species a has a conventional form, 



X [faip) Mk) [l±/c(p')] [l±/<i(fc')] - /c(p') /d(fc') [l±/a(p)] [l±/6(fe)]} • (2.2) 

We have introduced Us as the number of spin times color states for species s. {So Us equals 
6 for each quark or antiquark, and 16 for gluons.) Capital letters denote 4-vectors. The 
on-shell 4-momenta appearing inside the delta-function are to be understood as null vectors, 

P° = \p\ , etc. (2.3) 

We are using to denote Lorentz invariant momentum integration. 



ab 

cd 



The first term in curly braces in (2.2) is a loss term, and the second is a gain term. Ai 
denotes an effective scattering amplitude for the process ab cd, defined with a relativistic 
normalization for single particle states; its square, jA^cdP' should be understood as summed, 
not averaged, over the spins and colors of all four excitations (hence the prefactor of 1/ Ua) . 
The initial factor of 1/ (4|p|) is a combination of a final (or initial) state symmetry factor^^ of 
I together with the l/(2|p|) from the relativistic normalization of the scattering amplitude. 
Symmetry under time-reversal and particle interchange imply that 

iMt^ip, fc; p', fc') r = \Mtip, k; k', p') \^ = \M^,{k, p- p', k')\^ ^ \Mt,{p\ k'; p, k) \' , (2.5) 

etc. The effective scattering amplitude will itself be a functional of the distribution functions, 
since the density of other particles in the plasma determines screening lengths which affect 
the amplitude for soft scattering. This will be discussed explicitly in the next section. 

As Eq. (2.3) makes explicit, we have neglected medium-dependent corrections to quasipar- 
ticle dispersion relations in the overall kinematics of the collision terms, and in the particle 
velocity appearing in the convective derivative on the left side of the Boltzmann equation 
(2.1). Given our assumed separation of scales, these are relative 0{g^°') perturbations to the 
energy or velocity of a hard quasiparticle. Because we have assumed that distribution func- 
tions and observables are smooth functions on phase space, including (or excluding) these 
medium-dependent dispersion relation corrections will only affect subleading corrections to 
observables of interest. (This is discussed further in sections III and VI.) 

Now consider "1 2" processes. If isolated 1^2 processes are kinematically allowed by 
the effective thermal masses of the particles involved, and if there were no need to consider 
1 + N -ir^ 2 + N processes, then the appropriate 1 <-> 2 collision term would have a form 



When the species c and d are identical, a symmetry factor is required. When they are distinct, the final 
state is double-counted in the sum J2cd species. Hence a factor of ^ is needed in either case. 
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completely analogous to the 2 <-> 2 collision term: 

X [fM[l ± ± fc{k')] - h{p')fc{k')[l ± 




X {/a(p)A(fc)[l ± /e(p')] - /c(p')[l ± /a(p)][l ± ^(fc)]} ■ (2-6) 

With strictly massless kinematics, it is impossible to satisfy both energy and momentum 
conservation in a 1 2 particle process unless all particles are exactly coUinear. For small 
masses (compared to the energies of the primaries), the particles will be very nearly coUinear. 
One could then integrate over the small transverse momenta associated with the splitting 
to recast the collision term in the form 

c2-'^' [/] = E r -P'- 7." k'p) 

^\P\ 7o 

X ± Mj^'pW ± Uk'p)] - Mp'p)fc{k'p)[l ± fa{p)]} 

+ E / "^p' ^(\p\ + k-p') Yab{p'p;p, kp) 

X {fa{p)Mkp)[l±fc{p'p)]-fc{p'p)[l±fa{pm±Mkp)]}, (2.7) 

where we have ignored the small deviations from exact coUinearity when evaluating the dis- 
tribution functions. The factor of 7bc(p; p'p, kp) in the integrand is simply a way of param- 
eterizing the differential rate dT / dp dp' dk dflp for an a ^ 5c splitting processes, integrated 
over transverse momenta, excluding distribution functions and the longitudinal-momentum 
conserving (5-function. 

Any nearly-coUinear 2" process, now including 1 -|- A'^ 2 -|- soft scattering 

with emission events, can be cast into the general form of (2.7). The only difference is 
that the differential splitting/joining rates will now implicitly depend on the distribution 
functions for the A^ scatterers. All of the phase space integrations for those scatterers, plus 
the summation over A^, will be packaged into 7^^. The appropriate values of the splitting 
rates 7^^ will be determined simply by requiring that the collision term (2.7) reproduce 
previous results in the hterature for the rates of "1 <-> 2" processes. 

In the collision term (2.7), we have written the momenta of the splitting (or joining) 
particles as though they were exactly coUinear and as though their energy was exactly 
conserved. These particles actually receive 0{mes) kicks to their momentum and energy 
due to the soft interactions with the other A^ particles participating in the near-coUinear 
1-l-A^ <-> 2-1- A" process. As mentioned earher, this leads to a separation in the directions of the 
splitting particles by angles of order rrics / Pprimiay Treating this process as a strictly coUinear 
1^2 body process, and neglecting the soft O^nics) momentum transfers to other particles 
in the system, is an acceptable approximation because of our assumption that distribution 
functions (and observables) are smooth functions of momenta. If this assumption were 



15 




FIG. 5: Lowest-order diagrams for all 2 <-> 2 particle scattering processes in a gauge theory with 
fermions. Solid lines denote fermions and wiggly lines are gauge bosons. Time may be regarded as 
running horizontally, either way, so a diagram such as (D) represents both qq — ^ gg and gg qq. 

not satisfied, then we could not factorize the collision term into a product of distribution 
functions and an effective transition rate, as done above. Given our assumed separation 
of scales, the relative error introduced by this idealization is at most 0{g°')^ and hence 
irrelevant in a leading-order treatment. The differential rates are to be understood as 
summed over spins and colors of all three participants. Their explicit form will be discussed 
in section V. 



III. 2^2 PARTICLE MATRIX ELEMENTS 



Tree-level diagrams for all 2 2 particle processes in a QCD-like theory arc shown in Fig- 
ure 5. Evaluating these diagrams in vacuum {i.e., neglecting medium-dependent self-energy 
corrections), squaring the resulting amplitudes, and summing over spins and colors yields 
the matrix elements shown in Table 11.^^ Of course, vacuum matrix elements do not cor- 
rectly describe the scattering of quasiparticles propagating through a medium. In principle, 
one should recalculate the diagrams of Fig. 5 including appropriate medium-dependent self- 
energy and vertex corrections. But with generic hard momenta, for which all Mandclstam 
variables are of comparable size, these corrections are 0{g^°') effects and hence ignorable in 
a leading-order treatment. Medium-dependent effects can give 0(1) corrections to matrix 
elements if any Mandelstam variable is 0{mlg). Such momentum regions are phase space 
suppressed by at least g^" relative to generic hard momenta. Consequently, momentum 
regions where s ~ — t ~ —u ~ m^g (corresponding to either soft or collinear incoming parti- 
cles) give parametrically suppressed contributions, and medium-dependent corrections can 
be ignored in these regions. This implies that medium-dependent corrections do not have 
to be included in terms with denominators of s^.^^ 



Entries in Table II agree with well known QED results and the SU(3) results of Combridge et. al. [14]. 
Spacetime dependence of distribution functions will not be indicated explicitly, but it is important to keep 
in mind our basic assumption that, in the vicinity of any point x, distributions have negligible spacetime 
variation on a scale of tform- All results in this and later sections are local; they apply within a region of 
size tform about any particular point x of interest, using distribution functions evaluated at the point x. 
If one were to include self-energy corrections in these ,s-channel contributions, then there is a poten- 
tial subtlety if these terms can go on-shell. This is discussed in section VIA, but the conclusion that 
medium-dependent corrections can be ignored in s-channel processes (when external particles are treated 
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ab cd 

qiq2 ^ qiq2 , 
qiq2 ^ qiq2 , 
qiq2 ^ qiq2 , 

qiq2 ^ qiq2 

qiqi ^ qiqi , 
qiqi qiqi 

qiqi ^ qiqi 
qiqi ^ 92^2 

qiqi ^99 

qi9 ^qig, 
qi9 ^ qi9 

99^99 



Mtfl9^ 



did (s^ + u" 



t2 



dA \ tl 



dA \ £ 
did f t' + u' 

8 dp Cp (7 + — I — 8 dp Cp Ca j 2 



+ 16dpCp Cp- 



+ 















2 y 


1- 

' st 



t u 



-8dp ( - + - I + 8cZp Cp Ca 
' s u 



+ 

t2 



16dACi(3-^-^-^ 



TABLE II: Squares of vacuum matrix elements for 2^2 particle processes in a QCD-like theory, 
summed over spins and colors of all four particles. Here qi and q2 represent quarks of distinct flavors, 
gi and q2 are the associated antiquarks, and g represents a gauge boson, dp and dA denote the 
dimensions of the fundamental and adjoint representations, respectively, while Cp and Ca are the 
corresponding quadratic Casimirs. In an SU(A^) theory with fundamental representation fcrmions, 
dp = Ca = N, Cp = {N'^-l)/{2N), and dA = For a U(l) theory, dp = ^a = Cp = 1 

and Ca = 0. For SU(2), dp = Ca = 2, Cp = 3/4, and dA = 3, while for SU(3), dp = Ca = 3, 
Cp = 4/3, and dA = 8. Terms with underlined denominators are sufficiently infrared sensitive that 
medium-dependent self-energy corrections must be included, as discussed in the text. 



More problematic are the regions where —t or —u are O(m^g) while s is large, corre- 
sponding to a small angle scattering between hard particles. In this case, some of the 
matrix elements in the table arc enhanced by 0{s/mlff) or 0(s^/m^g). In Table II, terms 
with singly-underlined denominators indicate such infrared-sensitive contributions arising 
from soft gluon exchange, while terms with double-underlined denominators indicate IR 
sensitive contributions from a soft exchanged fermion. It is these underlined terms in which 
medium-dependent effects must be incorporated. ^° 



as massless) is unchanged. 

Sharp-eyed readers will notice that the s'^/{tu) and u'^/{st) terms appearing in the qq ^ qq and qq ^ qq 
matrix elements (squared) are not underlined, despite the fact that these contributions are infrared- 
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For the soft gluon exchange terms one must, in effect, reevaluate the small t (or small 
u) region of diagrams (A), (B), or (C) with the free gluon propagator on the internal line 
replaced by the appropriate non-equilibrium retarded gluon propagator, 

= — CHQ) - [Q' + nRet(g)];; . (3.1) 

(We have chosen Feynman gauge for convenience, but this is not required.) The required 
retarded gluon self-energy Il'^^{Q) is discussed in the next section. Evaluating the propa- 
gator (3.1) requires, in general, a non-trivial matrix inversion. (The matrix inversion may 
be performed explicitly in the special case of isotropic distributions; see Appendix A.) 

Because the self-energy only matters when the exchange momentum Q is soft, one has con- 
siderable freedom in precisely how the substitution (3.1) is implemented. Different choices, 
all equally acceptable for a leading-order treatment, include the following: 

1. Fully reevaluate the gluon exchange diagrams (A), (B), and (C), using the non- 
equilibrium propagator (3.1) for the internal gluon line. 

2. Introduce a separation scale /i satisfying me« -C -C Phard, and replace the free gluon 
propagator by the corrected propagator (3.1) only when Q"^ < [x^. 

3. Exploit the fact that soft gluon exchange between hard particles is spin-independent 
(to leading order) [8] . Write the IR sensitive matrix elements as the result one would 
have with fictitious scalar quarks, plus an IR insensitive remainder. Replace the IR 
sensitive part by the correct result for scalar quarks with medium corrections included. 

The final choice is technically the most convenient. It simply amounts to writing 

{s'W)/e = I + i {s-uf/e , su/e = I - 1 {s-uf/e , (3.2) 

and then replacing 



[s-u? 



\G^'\P-P')^,{P+PY{K+K'y\ ■ (3.3) 



To understand this, note that the square of the vacuum amplitude for i-channel gluon 
exchange between massless scalars is 



1^(0) (p-p%, (p+py {K+K'Y 



{P+P') ■ {K+K') 



{p-p'Y 



(3.4) 



For the \/v? terms, apply the same procedure with P' K' , which interchanges t and u. 



sensitive, albeit less so than the or terms. The s^/{tu) and u^/{st) terms arise from interference 
between t channel and either u or s channel gluon exchanges. They are sufficiently infrared singular to 
generate logarithmically divergent rates in the individual gain and loss parts of the collision term, but 
this log divergence (which would be cutoff by medium effects) cancels when the gain and loss rates are 
combined. Ifence, the apparent IR sensitivity of these terms may be ignored. 
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For the soft fermion exchange terms, one must similarly reevaluate the small t (or u) 
region of diagrams (D) and (E) with the internal free fermion propagator replaced by the 
non-equilibrium retarded fermion propagator, 

^ = -^^[(^-MQr'--§, (3.5) 

where = Q'^ — T,'^^^{Q). The retarded fermion self-energy ^Ret(Q) = Ifi^RetiQ) 
discussed in the next section. In the qq ^ gg matrix element, the net effect is to replace 

u 4Re[{P-Q){K-Qy]+sQ-Q* 

with — P^^—P'^^ — Ep^g^(P— P'), along with the corresponding replacement with P' <-> K' 
for the tu/v? term. For the qg <-> qg matrix element, the analogous replacement is 

s -AM{P-Q){P'-Qy]+tQ-Q* 
with Q^' now equaling pi^-K'^ - J:'^^{P-K'). Note that, in general, |Q • Q| V (Q • Q*)^- 



IV. SOFT SCREENING AND HARD EFFECTIVE MASSES 



To describe correctly the screening effects which cut off long range Coulomb interactions, 
we need the non-equilibrium generalization of the standard hard thermal loop (HTL) result 
for the retarded gauge boson self-energy n^etl^) with soft momentum, Q = O(meff). The 
general result for the non-equilibrium case has been previously derived by Mrowczynski and 
Thoma [15], who obtain^^ 



KM) 



2|p|(27r)3 dP^ 



-P^g^" + 



P-Q-ie 



P°=\p\ 



(4.1) 



where £ is a positive infinitesimal.^^ In this expression, the derivative df{p)/dP^ should 
be understood as zero. The sum runs over all species of excitations {i.e., g, u, u, d, d, ...), 
c^A is the dimension of the adjoint representation, and Cs denotes the quadratic Casimir 
in the color representation appropriate for species s.^^ (See the caption of Table II for 
specializations.) The self-energy (4.1) can also be rewritten in the form^'^ 

fsm 9' 



2\p\ (27r)3 



P-Q-ie 



(4.2) 



P°=\p\ 



The normalization of our distribution functions differ from Ref. [15] by a factor of 2, which appears in 
our expression as the spin degeneracy included in the factor z/g. To apply (4.1) for momenta Q of order 
meff, there is an implicit requirement that the background distributions f{x,p,t) not vary significantly 
on time or distance scales of order l/mcs. This is implied by our basic assumption in this paper that 
distributions are smooth on the still longer scale of tform- 
We use a ( — metric. 

For Abelian theories, replace g'^ Cs by the charge (squared) of the corresponding particle. 
To obtain this, replace (Pp/2\p\ in expression (4.1) by d^P 5{P^) 8(P"). Then integrate by parts, noting 
that the contribution where the derivative hits the delta function generates a factor of P'^, which vanishes 
when combined with the bracketed expression in (4.1). Then perform the integral to get back to d^p. 
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which is manifestly symmetric in the indices // and v. Simpler expressions may be given 
in the special case of isotropic distributions [16], where nRet(Q) has the same form as the 
equilibrium HTL result, but with a value of the Debye mass proportional to the integral J 
[Eq. (1.6)] over the distribution This is summarized in Appendix A. 

We mention in passing that the result (4.1) for nRet(Q) makes the implicit assumption 
that distributions f{p) do not vary significantly for changes of momentum of order q [which 
is Oirrics) in the case of interest]. Specifically, the derivation of (4.1) assumes that f{p+q) — 
f{p) can be approximated as q ■ 'Vpf{p). 

We also need the corresponding self-energy SRct(<5) for fermions with soft momentum Q, 
because this cuts off the small momentum-transfer behavior of fermion exchange diagrams 
hke Fig. 2. This was also derived by Mrowczynski and Thoma [15] but their result has a 
factor of 4 error. The corrected expression for a fermion of flavor s is ^■Ret,s{Q) — In ^Ret,s(Q) 
with 



P-Q-ie 



(4.3) 

PO = \p\ 



where /g, fs, and fs are the gluon, fermion, and anti-fermion distributions respectively (as 
always, per helicity and color state). The simple gamma- matrix structure ^ is a consequence 
of the chiral symmetry which results from the neglect of zero temperature fermion masses 
(relative to their medium-dependent effective mass). 

In addition to these soft self-energies, we will also need medium-dependent corrections 
to dispersion relations for hard particles that arc nearly on-shell. These will enter in the 
formulas that determine the near-coUinear "1 2" rates. As noted in the introduction, the 
resulting dispersion relation corrections for hard excitations, at leading order, turn out to 
take the simple form 

Q' + mls^O. (4.4) 

with a medium-dependent mass mcs- An efficient way to derive the values of these masses 
is by using the previous results for soft self- energies, which are valid for Q <^ Pscrccm in the 
intermediate regime nics <^ Q <^ Pscreen, where both the HTL approximation and the hard 
dispersion relation (4.4) are valid. 

Consider hard gauge bosons first. Only transverse polarizations are relevant because the 
longitudinal polarization decouples for hard momenta. We can use the massless dispersion 
relation = when evaluating the self-energy (4.2) for the purpose of obtaining the 
first correction to the dispersion relation. Working in light-cone coordinates defined by the 
direction of q, the second term in brackets in (4.2) is then proportional to 



Some readers may wonder whether there can be 0((5/pscreen) effects in the self-energy that were dropped 
in the HTL approximation but which affect the hard dispersion relation for Q > Pscreen, and which would 
cause meff to be a non- trivial function of Q/pscreen instead of a constant. One can check by explicit 
diagrammatic analysis of self-energies that, for = 0, this does not happen for fermions, scalars, or 
transverse gauge bosons. 
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For transverse choices of ji and the derivative vanishes. So the effective mass mefr,g of 
hard (transverse) gauge bosons comes only from the first term of (4.2), giving 
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I 



S fsip) . 



(4.6) 



s 



2\p\ (27r) 



Now consider hard fermions 



The effective Dirac equation is 



[<^-?Rct(g)]V' = o. 



(4.7) 



Multiplying on the left by another factor of — ^Ret) gives the condition 



{Q — SRet)^ — — 2Q ■ Ejiet + '^Kat — i 



(4.8) 



so that, to leading non-trivial order, the on-shell dispersion relation for hard fermions is 



From the result (4.3) for the self-energy, we immediately find that the effective mass for a 
hard fermion of flavor s is given by 



Note that the hard effective masses (4.6) and (4.10) are independent of the direction 
of the momentum q of the excitation, as well as its spin, even in the presence of general 
anisotropic distributions f{p). This is in contrast to soft screening, since the sclf-cncrgy 
n^g^(Q), given by (4.1), generically depends on the direction of q if distributions are not 
isotropic. The spin independence and isotropy of these (leading order) effective masses holds 
provided only that distribution functions are not themselves polarized. [This assumption 
underlies the HTL results (4.1) and (4.3).] This is special feature of hot gauge theories; 
in a generic anisotropic medium, no symmetry argument prevents splitting of dispersion 
relations into different branches depending on the spin of an excitation. In other words, a 
generic anisotropic medium is birefringent. 

If the self-energies (4.1) or (4.3) had led to spin-dependent dispersion relations for hard 
excitations, then even if distributions were not polarized at some initial time, the subsequent 
evolution of excitations through the birefringent medium would generate spin asymmetries. 
For such a system, it would be inconsistent to formulate an effective theory with spin in- 
dependent but anisotropic distribution functions. Fortunately, hot gauge theories do not 
behave this way. 

V. "1 ^ 2" PARTICLE PROCESSES 
A. Basic formulas 

The appropriate splitting/joining rates 7^^ characterizing near-collinear "1 ^ 2" processes 
may be extracted from Ref. [8]. Let n be a unit vector in the direction of propagation of 



g2 = 2g-ERet(<5)- 



(4.9) 




(4.10) 
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the splitting (or merging) hard particles, so that p — ph, p' — p' n, and k — kh. Then the 
required color and spin-summed effective matrix elements, consistently incorporating the 
LPM effect at leading order, may be expressed as 



p'2 _|_ p2 



lqgipn;p'n, kh) = -flg{ph;p'h, kn) = ^T^p-p ^q{P,P', k) , (5.1a) 

p'^+P^ + k^ 
p'^ p3 ^3 g 



lUpn]p'n,kh) = p^^^q"(^,-p',P), (5.1b) 

/4 _|_ 4 I 

lggipn;p'n,kn) = ^ip,P',k), (5.1c) 



where 



^ip',P,k) ^ 1^ / ^ 2h.ReFrih;p',p,k) (5.2) 

and a = g^/^An). The function Fp{h;p' ,p, k), for fixed given values of p', p, k and n, 
depends on a two-dimensional vector h which is perpendicular to n. Fp is the solution to 
the linear integral equation 

2h^iSE{h; p', p, k) Frih; p',p,k) + J 27r 5{vf. ■ Q) < < {A,{Q) [A^iQ]]*} 

X |(C,-iCA) [Fr{h;p',p,k) - Fr{h-kqr,p',P,k)] 
+ ICa [F^{h-p\p,k) - Fr{h+p'qr,p',P,k)] 
+ [Ff (/i;p',p,^) - F^{h-pqr,p',P,k)] ] , (5.3) 



which sums up multiple interactions. The four-vector Vfi = {^,n) is a null vector in the 
direction of n, and the vector q± is the part of q perpendicular to n. Once again, dg and 
Cg are the dimension and quadratic Casimir of the color representation for species s. The 
energy difference 5E is defined as 

SE[h;p.p.k) = — + — -— + — (6.4) 

and represents the energy denominator eg{k) + es{p) — ^sip') which appears in a p' pk 
splitting process. The variable h is related to the transverse momentum; see Ref. [8] for 
details. We will discuss momentarily the required correlator ((A^((5)[^i/((5)]*)) of the soft 
gauge field. 

Ref. [8], from which the above formulas for 7^^ were extracted, culminated in the derivation 
of the near-collinear "1 ^ 2" contribution to the total differential gluon production rate 
dTg/d^k for hard gluons with momentum k. In order to facilitate comparison with that 
reference, we show in Appendix B how to express the near-collinear contribution to dV^/d^k 
in terms of the 7^^ differential rates. 

The final element we need is the mean square fluctuations in soft momentum components 
of the gauge field in the medium, (l^An{Q)[Ai,{Q)]*J) . Formally, this is the Fourier transform 
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of the (non-equilibrium) HTL approximation to the Wightman gauge field correlator, ex- 
cluding the momentum and color conservation delta functions, 

I <fx d'y e'W'-^-Q-) (A^ix) At{y)} ^ {2n)'d^'\Q-Q') S^' {A,{Q)[A^{Q)r} . (5.5) 

This correlator characterizes the stochastic background fluctuations in which coUinear split- 
ting processes take place (see Ref. [8]). We are interested in the correlator for space-hke 
4- momenta Q, and physically (at leading order) it represents the correlation of the screened 
color fields carried by on-shell hard particles streaming randomly through the plasma. Hence, 
it may be understood as the (absolute) square of the amplitude shown in Fig. 6 integrated 
over the phase space of the hard particle, where the gauge propagator is to be understood 
as including the medium-dependent self-energy. This leads to 

{{A,mMQ)Y)) = G^^iQ) K^iQ) [G^^\Q)Y, (5.6) 
where G^^*'{Q) is the retarded propagator on the right side of (3.1) and 

K^{Q) = E V y '^^'^^ ^'^^^ ^ ^^^^^^^^ " " ^^-^^ 

s 

for soft momenta Q. Here, Vk = (1, k) is the 4-velocity of a hard particle with momentum k 
and g gives the current of this particle up to group factors. The /'s appear as initial state 
distributions and final-state Bose enhancement or Fermi blocking factors. For soft momenta 
Q, expression (5.7) may be simplified to give 

ni2"(g)=E'^^V J (^<^^^^^^)tl=^^^^^)]27r5(^fc-g). (5.8) 

One may alternatively derive this by summing hard thermal loops into the propagator of the 
Schwinger-Keldysh formalism, one of whose components gives the Wightman correlator. 
This is the origin of our notation H12 above, which is the one-loop Wightman current-current 
correlator and is an off-diagonal component of the self-energy in this formalism. 

In the special case of isotropic distribution functions, one can simplify substantially the 
expression (5.8) for Ili2{Q). Moreover, one can analytically reduce the d^Q integral appear- 
ing in the integral equation (5.3) to a two-dimensional integral over q±. See Appendix A 
for details. 



Not to be confused with time-ordered or retarded correlators of the gauge field. Once again, our basic 
assumption is that distribution functions are smooth on a time and distance scale of iform associated with 
the duration of a near-collinear "1 2" process. Hence, for the purpose of evaluating the non-equilibrium 
correlator in (5.5) for momenta of order nies, distribution functions may be treated as a;-independent. 
Spacetime variation in the non-equilibrium state will, of course, smear out the momentum conserving 
delta function in (5.5), but only by an amount which is irrelevant for our leading order treatment. 
For a general introduction, see chapter X of Ref. [17]. See also, for example, Eqs. (A34) and (A23) of 
Ref. [18]. 
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FIG. 6: Amplitude for the gauge field created by an on-shell hard particle in the medium. 

In equilibrium, the Wightman correlator is related to the retarded correlator by the 
fluctuation-dissipation theorem, which gives 

((^^(g)[A,(g)]*))equmbrium = 2 [n(g°) + 1] Im G)^,^* (Q) 

= -2[n(g°) + 1] G^^\Q) [lmUt,m [G^piQ)]*, (5-9) 

where n(e) = [e'^^— l]"-'^ is the usual equilibrium Bose distribution. For soft Q, the prefactor 

2[n{q^) + 1] can be replaced by 2T/q^. It is instructive to compare the non-equilibrium 
formula (5.8) for U^iQ) with that for ImU^etiQ)- The latter, for soft Q, can be extracted 
from the earlier expression (4.1) which gives 

^^<^iQ) = Y.''^^-df J J^3^"yki<l-^Mk)]n6{v,-Q). (5.10) 

The fluctuation-dissipation theorem (5.9) is satisfled in equilibrium by the results (5.8) and 
(5.10) for soft Q because equilibrium Bose or Fermi distribution functions ^^(e) = [e'^^^l]"^ 
satisfy 

q ■ Vfen(e) ^-fivq n,{e) [1 ± n,(e)] . (5.11) 
B. The formation time 

Earlier, we promised an explanation of the parametric formulas (1.10) and (1.11) for the 
formation time tform of a near-collinear "1 ^ 2" process and for the number iVform of soft 
collisions that take place in that time. We will give here a brief, superficial review, using the 
notation wc have adopted in this paper. For a more thorough discussion of the scales set 
by the LPM effect, the reader should consult the literature, such as Refs. [19-21]. We could 
deduce the scales by discussing the qualitative behavior of solutions of the actual equations 
(5.2) and (5.3) which incorporate the required physics, much as we did in the context of 
photo-emission in Ref. [11]. Instead, however, we will give here a more physical discussion 
which leads to the same results. 

One way to understand the basic scales is to begin by considering classical (soft) 
bremsstrahlung radiation of photons (rather than gluons), with wave number k^, emitted by 
a classical charged particle moving very close to the speed of light that undergoes random 
small-angle collisions. (This was the original classical picture of Landau and Pomeranchuk 
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[22].) Let T be the mean time between those coUisions and 9i the typical angle of deflection 
from each collision. The total deflection is then 9 ~ ^/N 9i, which we shall assume is small. 
If r is very large, there will on average be no interference between the bremsstrahlung fields 
created by successive collisions (for fixed k^). If r is very small, the bremsstrahlung field 
will not be able to resolve the individual collisions, and the field will be the same as that 
from a single coUision by the total angle 9. Since bremsstrahlung is at most logarithmically 
sensitive to the scattering angle 9, this smearing of N collisions into one collision will reduce 
the power radiated at this wavenumber by a factor of roughly N compared to what it would 
be if each collision could be treated independently. That is the LPM effect. 

The classical bremsstrahlung field from a scattering by angle 9 is dominated by radiation 
inside a cone of angle roughly 9. The first and last scatterings, at space-time points Xi 
and xn, will interfere significantly together if the phases in the factors exp(i^^^ ■ xi) and 
exp(i^^^ ■ X2) are comparable. For random collisions, that will happen if ■ {xn—xi) -C 1, 
which, for small 9, gives^^ 

A;^A^t(1 -cos^) ~ /c^A^T^^ < 1. (5.12) 

We have taken x%—Xi ~ |cciv— ~ Nr since the particle moves on nearly a straight-line 
trajectory at the speed of hght. For kjNT9'^ ^ 1, in contrast, the bremsstrahlung produced 
by the first and last collisions will be independent. The crossover criterion k^NT9'^ ~ 1, with 
9^ ~ N9f, then determines the typical number of collisions which are effectively smeared 
together by the LPM effect: 

N 1 1 ^ k-yT9l , (5.13) 

except that N must always be at least 1, since there must be a scattering to produce classical 
bremsstrahlung. 

A classical treatment of radiation breaks down when A;^ is no longer small compared to 
the energy E of the radiating charged particle. However, parametric estimates (as opposed 
to precise classical formulas) are still valid where the classical treatment first begins to break 
down. So we can still use the estimate (5.13) for N when k^ E, provided E — kj is not 
parametrically small compared to E. In a single soft collision of a hard particle with energy 
E, the typical defiection angle is 

(5.14) 

where q± is the typical transverse momentum transfer. Inserting this into (5.13) and setting 
kj E gives 

N ~ max 1^1, ^jE/irql)^ . (5.15) 

This estimate for hard bremsstrahlung applies equally well to the case of gluon emission. 
Setting q± ~ meff gives the estimate (1-11) previously quoted for A^form- 



For sharp single collisions (A''=l), there are actually two angular scales: the deflection angle 6 of the 

charged particle and the angle 6-^ that kj makes with the initial or final directions of the particle (whichever 
is smaller). There can then be logarithms arising from considering -C ^ in the bremsstrahlung rate, 
which we ignore. 
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^1 ~ 77^ ~ 77^ ^ 72- ' (5-16) 



When the formation time is large compared to r, then it is simply Nt by the definitions of 
and r. But if r ^ E/q\^ so that emission from different soft collisions do not significantly 
interfere, then the formation time is the time scale ti associated with bremsstrahlung from 
a single isolated collision. Classically, this is the time Ax^ for which ■ Ax ~ 1, which 
gives k^tiOf ~ 1 and 

E_ 

for /c-y ~ E. This same result can be found by examining how far off-shell in energy the 
internal hard line is in the basic bremsstrahlung processes of Fig. 3. One can nicely combine 
both cases in the single formula 

E 

tform ~ , (5.17) 

which, upon taking q± ~ rries, gives the earlier quoted result (1.10).^^ 

We specialized tok^ ^ E above. In this regime, the above estimates are equally applicable 
to photon or gluon emission. The behavior of the formation time for <^ E is not critical 
to understanding the conditions for applying our effective theory to the evolution of hard 
primaries. The dominant energy loss mechanism for primaries is via hard gluon emission 
processes with kg E rather than kg <^ E soft emissions. However, it is interesting to 
note that for soft emission, the case of photon emission is qualitatively different from gluon 
emission. Eqs. (5.13) and (5.16) show that the formation time for photons is much longer 
for k^ <^ E than it is for k^ ~ E. This is not true for gluon emission. Because the gluon 
can scatter by strong interactions, it cannot maintain its coherence over as long a time 
scale as a photon can. Soft scatterings involving the emitted gluon can change its direction 
by angles of order q±/kg ~ m^e/kg, which increase as kg decreases. Consequently, gluon 
bremsstrahlung with kg <^ E is actually less coherent than gluon bremsstrahlung with kg 
of order E/2. For gluon emission where the LPM effect is significant, the longest formation 
time is for the case where the energies of the two final particles are comparable. 



VI. DISCUSSION 

A. Dispersion relation corrections and double counting 

In section III we asserted that medium-dependent dispersion relation corrections on ex- 
ternal lines are, for hard particles, sub-leading corrections which may be neglected. There is. 



Replacing q± by in the estimates (5.15) and (5.17) amounts to an implicit assumption that the only 
important collisions are soft collisions with momentum transfers of order rrieff. Harder collisions are rarer 
than soft collisions. However, as noted in footnote 5, a single collision by an angle 6 ^ ^fN nies/E is no 
rarer than A'' consecutive soft collisions, each by angle nies/E. For large N, the same could be said of, for 
example, N/IQ consecutive collisions, each by angle ^/lOnies/E. This multiplicity of possibilities turns 
out to result in logarithmic corrections to the above analysis. Throughout this paper, we have consistently 
ignored logarithms in parametric estimates, and we continue to do so here. We have also ignored the effect 
of effective thermal masses, which cause hard particles to move slightly slower than the speed of light. 
Foi' Q-i- ^ mes, however, this does not affect any of our parametric estimates. 
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FIG. 7: The f-channel diagram for qg — > qg, drawn as a time-ordered diagram with time running 
from left to right. 

however, a potential subtlety concerning whether or not the internal line in a 2 2 process 
is kinematically allowed to go on-shell. If this occurs, then the 2 -^-^ 2 particle coUision rate 
will include the contribution from an on-shell 2 — > 1 process followed by a subsequent 1 — > 2 
process. Given the presence of explicit 1^2 particle collision terms in our effective theory, 
this would be inappropriate double-counting of the underlying scattering events. 

Consider, for example, the t-channcl qg qg process illustrated in Fig. 7. To incorporate 
a correct treatment of small angle scattering, as discussed in section III, one must include 
the HTL fermion self-energy (4.3) on the internal quark line. If medium-dependent effective 
masses are included on the external lines, then the internal quark line can go on-shell if 
mcfr,g > 2mefr,q. This condition is not satisfied in equilibrium QCD, but it can be satisfied 
with non-equilibrium distribiitions. It can also be satisfied, in equilibrium, in certain QCD- 
like theories. '^^ If this condition is satisfied, then the on-shell pole in the fermion propagator 
will generate a divergence in the two body coUision term C^^^. This divergence arises from 
an integration over the time difference between the creation and destruction of the virtual 
intermediate fermion, and refiects the fact that a calculation which just includes the HTL 
self-energy (4.3) is modeling that excitation as having an infinite lifetime. This divergence 
is, of course, unphysical and would effectively be replaced by the transport mean free time 
if further interactions with the medium were properly included. 

This situation is depicted in more detail (but still qualitatively) in Fig. 8, which shows the 
contribution to the loss part of the collision term C*^^^ (for hard particles) from t-channel 
qg ^ qg scattering, as a function of the invariant momentum transfer. For simplicity of 
presentation, this figure assumes that the momenta of scatterers, screeners and primaries 
are all within an 0(1) factor of a common scale T, and that the phase space distributions 
of all excitations on these scales are 0(1). In other words, the system is at most an 0(1) 
deviation from equilibrium. We have found it convenient to use r = t/^,/\^\ rather than t, 
and have sketched dCf^^ /dr vs. r. If the external particles were massless, then t (and hence 



For example, equilibrium SU(3) gauge theories with iVf Dirac fermions have TOeff,g = (^ + tk-^f)^^^ S*-^ 
and meff,q = gT. Hence meff,g > 2mcff,q if iVf > 10. 

This is because the HTL self-energies (4.1) and (4.3) have no imaginary part for timelike 4-momenta Q. 
One natural sounding but inadequate solution is to include the full thermal width on the intermediate line. 
However, this width is logarithmically IR divergent (in perturbation theory) due to sensitivity to long 
wavelength non-perturbative fluctuations in the gauge field. Moreover, merely including a width without 
simultaneously including additional interactions with the medium incorporates the wrong physics; the 
width is dominated by soft scattering, but a soft scattering event does not prevent the on-shell intermediate 
particle from propagating a large distance before breaking up. 
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FIG. 8: A qualitative sketch of the purely i-channel contribution of Fig. 7 for qg <-> qg to the 
loss part of the 2 2 collision term for a hard gluon or quark, plotted as dC^^^ /dr vs. r where 
T = t/ . Interferences with other channels are ignored. The labels mq and mg are short-hand for 
the hard effective quark and gluon masses mcff,q and mcff,g, and we have assumed mefr,g > 2mefr,q. 
The other scales listed above (1, g^, T, etc.) only denote parametric orders. 



r) would always be negative. The behavior between r ~ —T and r ~ —gT is dC/dr ~ g'^/r 
and is responsible for a leading-log contribution to the collision term. The (exaggerated) 
fall-off depicted for r <S —T reflects the decrease of initial state distributions for momenta 
large compared to T. The fall-off just above r ~ —gT is due to screening (that is, due 
to the inclusion of the HTL self-energy for the internal quark line). For r < — gT, the 
contributions to Cj^t^^ are dominated by hard initial particles whose trajectories intersect at 
an angle 612 of order 1. Each particle is deflected in the collision by an angle of order \t\/T. 
The contributions in this kinematic region are not sensitive, at leading order, to whether or 
not one uses massive or massless dispersion relations for the external lines. In contrast, for 
r > 0, the contributions are dominated by hard initial particles that are nearly collinear, 
with 612 ~ g. The peak at r = nics^q represents the nearly- collinear process of Fig. 7. 
[Kinematics forces this process to be nearly collinear because the 0{gT) masses mefr,q and 
mefT,g are small compared to the hard 0{T) momenta.] 

Prom the plot, one can see that two regions make signiflcant contributions to the collision 
term. The flrst region is — T < r < gT, which reflects genuine 2^2 processes and makes an 
0{g^T) contribution to Cioss- The second region is |t— meif,q| < g^T, which is double counting 
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FIG. 9: s-channel diagram for qq ^ qq. 




FIG. 10: A qualitative sketch of the purely s-channel contribution of Fig. 9 for qq ^ qq to 
the collision term C for a hard quark, plotted as dC/d'^/s vs. ^/s. The dot-dash line shows the 
dependence if meff,q and meff,g are both set to zero. 



and incorrectly treating the 1^2 processes that arc described by the "1 ^ 2" collision term. 
Fortunately, our treatment of "1 2" processes correctly handles off-shellness in energy 
as large as the inverse formation time 0{g'^T) (see section I A) and so already correctly 
accounts for the physics in this region where the energy is off-shell by only < g^T. Note 
that r's further out on the peak at nies^q in Fig. 1 (for instance, |r — mefr,q| ~ g^T), give 
a contribution to Cioss from this particular diagram that is subleading compared to genuine 
2^2 contributions and which may therefore be ignored. 

To formulate a correct collision term which only includes "genuine" 2 -^-^ 2 processes, 
one must somehow keep the contribution from the first region and eliminate the second. 
Although there are many ways one could accomplish this, the simplest solution is to include 
the HTL self-energy on internal lines (where it is needed to describe screening) but to treat 
external lines as massless. 

Exactly the same issue can arise with s-channel processes, as illustrated in Fig. 9 for 
qq <-> qq. A plot analogous to the one discussed above, but this time showing the contribution 
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to dC'^^^ /dy/s vs. y/s for this s-channel annihilation reaction, is shown in Fig. 10. If HTL 
self-energies are included in the internal gluon line, then the virtual gluon can go on-shell 
in this qq ^ g ^ qq process even when the external quarks are treated as massless, since 
two massless particles are kinematically allowed to combine to create one massive one. For 
s-channel processes, however, it is not necessary to include HTL self-energies on the internal 
line in the first place. Unlike the case of t and rt-channel processes, they are not required 
to control infrared sensitivity. So one solution which avoids all partial double counting of 
1^2 processes while not affecting the leading-order result for the true 2^2 contributions 
is to treat all external particles as massless, and include HTL self-energies in t-channel and 
w-channel internal propagators but not in s-channel ones. This is the approach presented in 
section HL The dot-dash hne hiding under the peak in Fig. 10 and smoothly decreasing at 
small s illustrates the result of this prescription for s-channel processes. 

There are alternative possibilities which are equally valid at leading order. For example, 
one could include HTL self-energies multiplied by the step function Q{Q'^) on all internal lines 
in 2 <-> 2 processes, so that they only affect spacelike propagators. With external particles 
treated as massless, this would also avoid double-counting mistakes. But attempting to 
"improve" the effective theory by including both medium-dependent self-energies on internal 
lines and dispersion corrections on external lines is simply wrong — unless the contributions 
from 2 2 processes degenerating into two independent scatterings are carefully separated 
and subtracted. Although this could be done consistently,^^ it is needlessly comphcated 
compared to the simple approach of treating external lines as massless and only inserting 
medium-dependent self-energies where they are truly required. 

B. Additional scattering processes 

One may wonder if any additional scattering processes need to be included in a leading- 
order effective theory. To examine this, first note that processes whose rates are paramet- 
rically slower than hard 2 2 or 1 2 processes will have negligible {i.e., subleading in 
g) effect on the dynamics of a quasiparticle. For example, consider adding an additional 
radiated gluon to a hard (i.e., large momentum transfer) 2 2 scattering process. The 
case of qq qqg is illustrated in Fig. IL In the context of high-energy collisions in vacuum, 
it is well-known that bremsstrahlung gluons cost a factor of g^ times logarithms associated 
with coUinear and soft infrared enhancements. Soft gluon emission is unimportant for our 
leading-order effective theory which only describes the dynamics of hard quasiparticles. And 
in any case, sensitivity to small gluon momenta will be cut off in a medium by the effective 
mass of the gluon. Collincar logarithms will similarly be cut off by the effective masses of 
the quark and gluon. However, in a medium there is also a [l+/g] final state statistical 
factor associated with the radiated gluon. This factor can be parametrically large. But, as 
discussed in section I C, our assumptions limit distribution functions to be small compared 
to 0{g~'^) for p ^ meflf. The upshot is that each radiated gluon suppresses the transition 



Exactly the same issue of 2 ^ 2 scattering processes degenerating into two independent 1^2 scatterings 
arises anytime one attempts to formulate a kinetic theory for unstable particles. Ref. [23], for example, 
includes a discussion of the need to subtract the degenerating part of 2 <-> 2 scattering rates in the context 
of various models of GUT-scale baryogenesis. 
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FIG. 11: A subleading hard qq qqg process. 

rate by a factor that is parametrically smalL [If we focus on bremsstrahlung of hard gluons, 
then condition (1.15) implies the suppression is at least g'^°', possibly times logs of 1/g.] 
Therefore, a primary quark or gluon will experience a parametrically large number of hard 
scatterings without gluon bremsstrahlung before it undergoes one with bremsstrahlung. 

The simple fact that a process is suppressed compared to others does not automatically 
make it irrelevant at leading order. For instance, 2^2 scattering does not change the 
total number of hard quasiparticles. Correctly including the dominant number changing 
processes, even if they are parametrically slow compared to the large angle scattering time, 
is necessary for a leading-order calculation of physics that depends on equilibration in the 
number of quasiparticles (such as bulk viscosity). However, bremsstrahlung from hard 
2 2 scattering is not the fastest number changing process. 

The most likely scattering events (as discussed in section I A) are soft scatterings with 
momentum transfer of order rries- The relative suppression for radiated bremsstrahlung 
gluons is the same as above (up to logarithms) and is at least g'^"' for hard bremsstrahlung 
gluons. That is, an excitation will experience a parametrically large number of soft scat- 
terings unaccompanied by hard gluon bremsstrahlung before it undergoes a soft scattering 
with bremsstrahlung. And soft scattering with single bremsstrahlung is the fastest number- 
changing process (illustrated in Fig. 3) for hard particles. Fortunately, this process has 
already been included in our effective kinetic theory. It is part of the A'^ -|- 1 — > A^ -|- 2 
processes which have been summed up in our effective "1 ^ 2" near-coUinear transition 
rates. 

The basic point is that a hard scattering accompanied by bremsstrahlung does not change 
the distributions of quarks or gluons in any way which is distinct from the much more 
rapid effects of hard 2 <-> 2 scatterings together with soft scattering accompanied by gluon 
bremsstrahlung. In contrast, a soft scattering accompanied by gluon bremsstrahlung cannot 
be neglected, both because it changes particle number and because a single scattering of this 
type can produce an 0(1) change in the momentum of an excitation (unlike the more rapid 
2 <-> 2 soft scatterings) at a rate which can be competitive with hard 2 <-> 2 scattering. 



This is discussed at some length, in the context of bulk viscosity for scalar theory, in Refs. [24, 25]. 
Recall that bremsstrahlung gluons are dominated by angles less than or order of the deflection angle in 
the underlying 2^2 scattering event. Diagrammatically, this occurs because of cancellations between 
different diagrams, such as those of Fig. 3. But it can also be seen by reviewing classical formulas for 
the intensity of bremsstrahlung radiation. Since the momentum transfer of the most prevalent collisions 
is order nies, the deflection angle is order meflf/pprimary for a primary excitation. In local equilibrium 
settings, this is order g. 
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FIG. 12: A diagram contributing to double gluon bremsstrahlung. 




FIG. 13: Double gluon bremsstrahlung with additional soft scattering. 

To make the above discussion more concrete, let us briefly specialize to typical excitations 
in near-equilibrium systems. In that case, the rate of hard scattering with hard gluon 
bremsstrahlung is 0{g^T), which is 0{g^) suppressed relative to the 0{g^T) rate of hard 
2 <-> 2 processes. The rate of soft scattering with hard bremsstrahlung is 0{g'^T), which is 
0{g'^) suppressed relative to the 0{g^T) rate of straight 2 <-> 2 soft scattering. 

As a further example of the suppression of higher-order processes, consider soft 2 2 
scattering with hard double bremsstrahlung, as depicted in Fig. 12, which we might call 
a "1 ^ 3" splitting. As already discussed, this would be suppressed by g"^" compared to 
single bremsstrahlung. It is therefore irrelevant at leading order, since a parametrically large 
number of "1 <-> 2" single bremsstrahlung events will occur for every one of these double 
bremsstrahlung events. [Near equilibrium, the rate for this "1 <-> 3" double bremsstrahlung 
is g^T, compared to the g^T rate of "1 ^ 2" processes.] 

Now consider adding one or more additional scatterings to the double bremsstrahlung 
process, to generate a diagram like the one shown in Fig. 13. Adding additional soft scat- 
terings can increase the cross-section, because the internal line running from x to |/ in the 
figure can go on shell.^^ This part of the amplitude, however, really represents two con- 
secutive "1 ^ 2" processes and so is already accounted for. The process of Fig. 13 differs 
significantly from two consecutive collisions only when the time between x and y is com- 
parable to the duration [l/g'^T near equilibrium) of the individual "1 ^ 2" processes. In 
momentum space, that corresponds to the kinematic region where that propagator is just 
as off-shell in energy as the intermediate quark line in the double bremsstrahlung, single 
soft scattering case of Fig. 12 considered previously. In other words, there is no additional 



For comparison, the internal quark lines in the double bremsstrahlung process of Fig. 12 have virtualities 
~ miff since the angle between each hard gluon and the emitting quark line is 0(meff/pprimary)- 
(See the previous footnote.) Hence these internal lines are off-shell in energy by an amount of order 
P'^/E m^jf/pprimary (or fif^T in the near-equilibrium case). These internal lines are prevented from 
being more on-shell by the medium-dependent effective masses of particles. 
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on-shell enhancement except when the process degenerates into separate scattering events, 
which are already included in the effective kinetic theory. 

In summary, multiple gluon bremsstrahlung processes are either indistinguishable from 
a sequence of "1 ^ 2" and 2^2 processes or else arc parametrically slower and do not 
accomplish any relevant relaxation not already provided by a sequence of faster processes. 
After examining these, and other possibilities, we are unaware of any processes which can 
affect reasonable obscrvablcs at leading order in g, beyond those which have already been 
included in our effective theory. 



C. Soft gauge field instabilities 

In section III, we noted that in processes involving soft gluon exchange, one must include 
the medium-dependent retarded self-energy in the intermediate gluon propagator in order to 
incorporate the screening (and Landau damping) of long range gauge field interactions. For 
situations involving substantial departures from equilibrium, there is an important question 
which has not yet been addressed — does the medium-dependent self-energy actually cut off 
long range interactions? In other words, are the resulting 2^2 non-equilibrium scattering 
rates well-defined? 

The relevant part of the phase space integral for soft gluon exchange processes involves 
an integral over the spacelike momentum transfer of the product of advanced and retarded 
non-equilibrium gluon propagators. 



/ 

Jq 



d'Q Q{Q') [Q' + nRet(g)]:,' [Q' + nRet(g)i;' . (e.i) 



This is then contracted with factors which, by gauge invariance, are necessarily transverse 
to Q. If x^iQ) aiid \a{Q) denote the eigenvectors and associated eigenvalues of Iljiei{Q), 
then the potentially dangerous part of this integral is the piece involving projections onto 
the same transverse eigenvector x"'{Q) from both propagators, 

.4^o.^2^ x'',{Q)x:{Qrxi{Q)x${Qr 

Q«p1 ^ ^ [Q' + Re A„(Q)]2 + [Im A„(Q)]^ " ^'"'^ 

The integrand will be singular if the gluon propagator has poles at real spacelike momenta, 
which can potentially happen if an eigenvalue Xa{Q) is real and negative for some domain 
of spacelike momenta. One may show that a spacelike pole in the gluon propagator creates 
a non-integrable singularity in the integral (6.2) and generates a logarithmic divergence in 
the soft collision rate. The presence of such a spacelike pole in the non-equilibrium retarded 
gluon propagator would imply that the corresponding modes of the soft gauge field have 



This assertion does deserve a few caveats. Since QCD exactly conserves the net fermion number in each 
flavor, but weak interactions do not, weak interaction collision terms cannot be neglected if one is interested 
in the evolution of flavor asymmetries on suflaciently long time scales. Similarly, in hot electroweak theory 

expUcit baryon prodTiction/dcstruction terms representing the effects of non-pcrturbativc baryon number 
changing processes can be relevant on time scales large compared to the mean free times of processes 
discussed in this paper. See Ref. [2] for more discussion of these points. 
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exponentially growing behavior in time. To our knoweledge, the first discussion of such 
instabilities in the context of QCD plasmas was by Mrowczyhski in Ref. [13]. 

In equilibrium, no such instability is possible. More generally, such instabilities do not 
appear if distribution functions are isotropic, but otherwise arbitrarily far from equilibrium. 
This is a consequence of the fact, discussed in Appendix A, that if distribution functions are 
rotationally invariant then the non-equilibrium HTL self-energies turn out to be proportional 
to their equilibrium form. 

For anisotropic but parity invariant distributions, it turns out that spacclikc poles are 
generically present. To see this, first note that the HTL self-energy (4.1) does not depend 
on the magnitude of Q, but only on its (four-dimensional) direction. And, as may be seen 
from Eq. (5.10), for any parity invariant distribution the imaginary part of IlRetiQ) vanishes 
identically when q'^ — 0, implying that the eigenvalues of the zero frequency self-energy 
are purely real. Within the = surface, if there are directions in which ReAa(Q) is 
negative then there will be singularities in the integral (6.2). For generic anisotropic but 
parity invariant distributions, there are such directions. One may show that the angular 
average of the trace of the HTL spatial self-energy, H-((5), vanishes at zero frequency, which 
means that the angular average of the sum of the two transverse eigenvalues of the spatial 
zero-frequency self-energy vanishes. Hence, if there is any direction in which an eigenvalue 
of the static spatial HTL self-energy is positive, then there must be some direction in which 
an eigenvalue is negative. The only way instabilities can be avoided is if the zero-frequency 
spatial gluon HTL self-energy vanishes identically, as it does for isotropic distributions. 

If one perturbs a parity invariant distribution by adding a parity non-invariant compo- 
nent, then a continuity argument shows that instabilities will still generically be present for 
sufficiently small deformations away from the parity invariant case. [In the presence of an 
arbitrarily small parity non-invariant perturbation, Im Xa{Q) must vanish on some surface 
which is a small deformation of the g° = plane. Within this surface, there must still be 
directions in which ReXa{Q) is negative, since there are such directions in the absence of the 
deformation.] Whether these spacelikc poles of the gluon propagator persist for completely 
general non-parity invariant, anisotropic distributions is not yet clear to us. 

For a given set of distribution functions, if spacclikc poles are present then the character- 
istic wave vector (and growth rate) of the associated soft gauge field instabilities will at most 
be of order of the non-equilibrium effective mass rries, since this is the scale which character- 
izes the size of medium-dependent self-energy corrections. For distribution functions which 
are parametrically close to isotropy, the wave vector and growth rate of soft instabilities will 
be parametrically smaller than rUes- If any gauge field modes with wave vectors of order 
ruefi are unstable, then the growth of such instabihties would be expected to lead to spatial 
(and temporal) inhomogeneities in distribution functions on length scales of order l/m^^. 
Consequently, the growth of such instabilities should lead to a violation of the spacetime 
smoothness condition which underlies our effective kinetic theory. The physics which cuts 
off the growth in these instabilities, and removes the divergence in the soft scattering rate, 
can only come from including the effects of spacetime inhomogeneities inside the evaluation 
of the effective scattering rate. This means one can no longer use the derivative expansion 
which underlies the effective kinetic theory. 

As far as our effective kinetic theory is concerned, the net result is that its domain of 
validity is smaller than anticipated. We initially required that anisotropics in distribution 
functions not be parametrically large. But it appears that 0(1) anisotropics in the excita- 
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tions responsible for screening will lead to violations of the assumed spacetime smoothness 
condition on a time scale that is shorter than the minimal time scale tform for which the ef- 
fective theory was designed. However, parametrically small anisotropics are still acceptable. 
An 0{g) anisotropy in the distribution of "screeners" can only generate instabilities whose 
wave vectors are parametrically small compared to rries- In a leading-order treatment, one 
may simply excise this "ultrasoft" momentum region from the phase space integral (6.2). 
The resulting ambiguity in the effective collision rates will be subleading in g, and our 
effective kinetic theory should remain valid for the intended class of observables — those 
primarily sensitive to the dynamics of hard excitations. 

Clearly, it would be interesting to study the effects of soft gauge field instabilities in non- 
equilibrium systems and try to understand their influence on physical observables which 
probe the relevant soft or very-soft dynamics. This is a topic for future work. 

D. Effective kinetic theory beyond leading order? 

Typical effective theories (such as heavy quark theory, or non-relativistic QED) can be 
systematically improved order-by- order in powers of the ratio of scales which underlies the 
effective theory. Having constructed a leading-order effective kinetic theory for the dynamics 
of a hot gauge theory, it is natural to ask whether one can formulate a beyond-leading-order 
kinetic theory which will correctly incorporate relative corrections suppressed by one or more 
powers of g. This is an interesting open question. 

Consider, for simplicity, the case of systems which differ from equilibrium by at most 
0(1), so that all relevant hard momenta are 0{T). Any attempt to construct an effective 
theory of hot dynamics beyond leading order must handle numerous different sources of 
subleading corrections. These include: 

1. Kinematic mass corrections of order ^eff/i^Lrd ^ 9^- These are everywhere: in the con- 
vective derivative of the Boltzmann equation, in the overall kinematics of the collision 
terms, inside the effective 2^2 scattering amplitudes, etc. Consistently including 
such corrections should be feasible, but will force one to separate and subtract the 
degenerating parts of scattering rates, as discussed in section VIA. 

2. Contributions from higher order tree processes, such as bremsstrahlung from hard 
scattering. Numerous such additional processes appear at 0{g'^) and would need to be 
included in the collision terms, again with appropriate care to eliminate phase space 
regions where an intermediate line goes on-shell and the process separates into multiple 
scattering events. 

3. Loop corrections to 2 ^ 2 effective scattering amplitudes. For hard scattering, these 
will be order g^ effects, but for soft scattering, the relevant loop expansion parameter 
is g, not g^. Also, the HTL approximation to the self-energies required in soft exchange 
processes receive 0{g) corrections because of kinematic approximations made in the 
HTL results. 

4. Subleading corrections to effective near-collinear transition rates, and a proper treat- 
ment of the soft emission region. We believe these enter at 0{g). Evaluating such 
corrections would require a next-to-leading order treatment of LPM suppression. This 
is unknown territory. 



35 



5. Contributions from soft {p ~ gT) on-shell excitations. The size of such contributions 
depends on the sensitivity of observables of interest to soft momenta. For observables 
Uke fermion current densities or the traceless part of the stress tensor (whose behaviors 
determines diffusion constants and shear viscosity), soft contributions are suppressed 
by 0{g'^) or more. Incorporating soft contributions requires formulating a kinetic 
theory which correctly describes both hard (ultrarelativistic) and soft (non-relativistic) 
excitations. 

6. Contributions from non-perturbative gauge field dynamics on the g'^T (ultrasoft) scale. 

The importance of very small angle scattering via ultrasoft gauge boson exchange is 
suppressed, relative to soft exchange, by g^; so we expect ultrasoft physics effects to 
enter as g^ corrections to our effective kinetic theory. This means that nonperturbative 
inputs will be necessary to formulate a kinetic theory which correctly describes 0{g^) 
corrections. We have no idea how this could be done in practice. 

7. Corrections due to the uncertainty in energy of excitations. The relative size of such 
corrections is controlled by the inverse of an excitation's energy times its mean free 
time between scatterings. For soft scatterings of hard excitations, this is 0{g^). Cor- 
rectly incorporating such quantum corrections to kinetic theory is an interesting open 
problem. 

Consistently incorporating 0{g) corrections may well be feasible, but extending the ef- 
fective kinetic theory to include 0{g'^) effects involves major conceptual challenges as well 
as technical difficulty. 

VII. CONCLUSION 

We have argued that quasiparticle dynamics in relativistic plasmas associated with hot, 
weakly-coupled gauge theories (such as QCD at asymptotically high temperature T) can be 
described by an effective kinetic theory, valid on sufficiently large time and distance scales. 
This effective theory is adequate for performing leading-order evaluations of observables 
(such as transport coefficients and energy loss rates) which are dominantly sensitive to the 
dynamics of typical ultrarelativistic excitations. In other words, our effective theory ne- 
glects effects which generate relative corrections suppressed by powers of the gauge coupling 
g, but correctly includes all orders in l/logg~^. To construct such a leading-order effective 
theory, it was necessary to include in the collision term of the kinetic theory both 2^2 
particle scattering processes as well as effective "1 <-> 2" coUinear splitting and merging 
processes which represent the net effect of nearly coUinear bremsstrahlung and pair produc- 
tion/annihilation processes taking place in the presence of fluctuations in the background 
gauge field. 

Our effective kinetic theory is applicable not only to near-equilibrium systems (relevant 
for the calculation of transport coefficients), but also to highly non-equilibrium situations, 
provided the distribution functions satisfy the conditions discussed in section I C [as amended 
in section VIC] which, in particular, require that there is a clear separation between the 
Debye screening scale and the momenta of typical excitations of interest, and that the 
excitations responsible for screening be close to isotropic. These conditions can be satisfied 
in asymptotically high temperature QCD, where the running coupling g{T) is truly small. 



36 



They may also be satisfied at intermediate stages of collisions between arbitrarily large 
nuclei at asymptotically high energies [3]. What, if any, utility this effective theory has for 
understanding real heavy ion collisions at accessible energies is not yet clear. However, we 
believe that understanding dynamics in weakly coupled asymptotic regimes is a necessary 
and useful prerequisite to understanding dynamics in more strongly coupled regimes. 
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APPENDIX A: ISOTROPIC DISTRIBUTIONS 



Isotropic (rotationally invariant) distributions represent a physically interesting class of 
problems intermediate between completely general non-equilibrium systems on one hand, 
and equilibrium or near-equilibrium systems on the other. It is widely appreciated that 
substantial simplifications occur for equilibrium systems. For isotropic systems, almost as 
high a level of simplification is possible. The following results are written in the plasma rest 
frame, which is uniquely defined for an isotropic system. 

We begin with the retarded fermion self-energy SRet(-^) = 7At^Ret(-^)' ^ 
momentum K = (k^,k) such that ~ Vi^ = A; <S Pscreen- With isotropic distributions, 
this soft self-energy has exactly the same structure as in equilibrium [26, 27], 
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(A2) 



The overall coefficient m^g f is given by the integral (4.10) which determines the effective 
mass for hard fermions. In the literature, this self-energy is more conventionally written in 



terms of Afq 



|mgg f, because Mq is then the frequency of oscillation for a fc = fermion. 



The case of the retarded gauge field self-energy nRet(-f^) is similar. As in thermal equi- 
librium [28, 29], the self-energy can be decomposed into longitudinal and transverse pieces. 



with 



fc2 ' 
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F [(A;0) 
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(A4) 
(A5) 



Here k^ = K'^ + u - K = (0, fc) denotes the part of the 4-momentum K orthogonal to the 
rest frame 4- velocity u. The projectors P^'^, Q^^ , and K^K'^/K'^ are mutually orthogonal 
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and sum to the metric, P'^" + Q'*'" + K^^K'^/K'^ — rj^". The transverse and longitudinal 
self-energies are 
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where ?7iggg is the asymptotic gluon mass defined in Eq. (4.6). Equivalent forms in the 
literature are more commonly written in terms of the leading order Debye mass = 2m^g ^ , 
or occasionally in terms of the leading order plasma frequency cjpj = m^/S. Note that the 
definition of IIl is not uniform in the literature (even in previous work by the authors of the 
present paper!). The above notation agrees with that of Weldon [29]. The other common 
usage is that of Braaten and Pisarski [30], who define IIl to be —k'^/K'^ times our value, so 

_j_ ■|--|-(Braaten-Pisarski) JX^^ 

Rotational symmetry, even in the absence of equilibrium, is sufficient to derive a relation 
between the imaginary part of the retarded self-energy and the Wightman self-energy, at 
soft momenta. Namely, 

9T 

where 
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and Xg and Jg are the integrals defined in Eqs. (1.7) and (1.6), evaluated with the distribution 
function for species s. Relation (AS) is just the small frequency form of the equilibrium 
fiuctuation-dissipation relation (5.9), but with the equilibrium temperature replaced by T*. 

Consequently, the Wightman correlation function, for soft momenta, is merely a rescaled 
version of the thermal Wightman correlation function. 



{{A,{K)[A^{K)r 
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at the temperature T^g for which the equilibrium Debye mass coincides with the correct 
effective Debye mass,^^ iTi^{Tes)\eqmi. — Sm^^g. Note that the value of T* is not in general 
the same as Tefj. 

Relation (AlO) permits one to use recent results of Aurenche et al. [31] to reduce the four 
dimensional integral involving the Wightman correlator appearing in the integral equation 
(5.3) down to a two dimensional integral over transverse momenta. One may show that 
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where h{q±) is any function of q±. 
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Explicitly, for an SU(-/Vc) theory with N{ fundamental Dirac fermions, ^{Nc+Nf Cp dp/dA) 9^T^s = '^eff.g- 
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APPENDIX B: RELATIONSHIP OF dT^/d^k TO 7^ 



a 



The leading- order equilibrium differential rate for production for hard gluons, as defined 
in Ref. [8], corresponds to the gain part of gluon collision terms (2.2) and (2.7), evaluated in 
equilibrium and multiplied by i/g/(27r)^. Explicitly, the near-coUinear LPM-suppressed part 
of the production rate for hard gluons with momentum k — kn is, at leading order. 



dr^™ _ 1 dTi 



d^k 47r/c2 dk 

= p U y dp dp' 5{p'+p-k) -f^gg{kn;p'h,ph) n^{p') nb(p) 

/•oo 

+ dp dp' 5{p'-p-k) -f^g{p'n;pn, kh) ny,{p') [l+nb(p)] 
^0 

roo 

+ N{ dp dp' S (p'+p-k) ^qq{kh; p'h, ph) nf{p') nf{p) 
Jo 

+ 2N{ dp dp' 6{p'-p-k) Ygg{p'n;pn, kn) nf(p') [l-rif (p)] j . (Bl) 

The i in the first term is an initial-state symmetry factor, the 2 multiplying the final term 
refiects the identical contributions from q ^ qg and q — > qg, and Nf is the number of Dirac 
fermion flavors. n]^{uj) and nf(a;) are equilibrium Bose and Fermi distribution functions, 
respectively. Ref. [8] expresses results more compactly by using crossing symmetries. For a 
more direct comparison with that paper, expression (Bl) can be rewritten as 

= il±^ dp 5{p'-p-k) ilg{p'h;ph, kh) n^ip') [1+nM 

+ Nf J dp dp' S{p'—p—k) j^g{p'ri;pri, kh) nf{p') [1— nf(p)]^ . (B2) 

However, the authors of Ref. [8] should be profoundly chastised for not pointing out 
that this differential gluon production rate is, in fact, an infrared divergent quantity. The 
problem arises from the p — > region of the g gg term in (B2). This portion of the 
integral represents processes in which a hard gluon with momentum p' nearly equal to k 
experiences a soft scattering with emission or absorption of a soft gluon to yield a hard gluon 
with momentum k. But physical quantities can only depend on this production rate minus 
the corresponding rate at which gluons are scattered out of the mode k, and the infrared 
sensitivity cancels in the difference of these rates. In other words, although the production 
rate dT^^^/d^k is not actually well-defined, the complete collision terms (2.7) built from the 
same near-coUinear transition amplitudes are infrared safe. 
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